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Abstract 



The parameters of inspiralling compact binaries can be estimated using 
matched filtering of gravitational-waveform templates against the output of 
laser-interferometric gravitational- wave detectors. The estimates are most 
sensitive to the accuracy with which the phases of the template and signal 
waveforms match over the many cycles received in the detector frequency 
bandwidth. Using a recently calculated formula, accurate to second post- 
Newtonian (2PN) order [order (t>/c) 4 , where v is the orbital velocity], for 
the frequency sweep (dF/dt) induced by gravitational radiation damping, we 
study the statistical errors in the determination of such source parameters as 
the "chirp mass" Ai, reduced mass //, and spin parameters (3 and a (related 
to spin-orbit and spin-spin effects, respectively). We find that previous results 
using template phasing accurate to 1.5PN order actually underestimated the 
errors in A4, fi, and 0. Templates with 2PN phasing yield somewhat larger 
measurement errors because the 2PN corrections act to suppress slightly the 
importance of spin-orbit contributions to the phase, thereby increasing the 
measurement error on (3. This, in turn, results in larger measurement errors 
on M and fi because of the strong correlations among the parameters. For 
two inspiralling neutron stars, the measurement errors increase by less than 
16 percent. 

PACS numbers: 04.25.Nx, 04.30.-w, 97.60. Jd, 97.60.Lf 
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I. INTRODUCTION 



Inspiralling compact binary systems, composed of neutron stars and/or black holes, have 
been identified [|I]|| as the most promising source of gravitational waves for interferometric 
detectors such as the American LIGO (Laser Interferometer Gravitational-wave Observatory 
@) and the European Virgo |4j]. These systems evolve under the influence of gravitational 
radiation reaction, so that the gravitational-wave signal increases in amplitude as its fre- 
quency sweeps through the detector frequency bandwidth, from approximately 10 Hz to 
1000 Hz. (This characteristic signal is often referred to as a "chirp.") Inspiralling compact 
binaries are especially promising because recent estimates [^,0 indicate that their event rate 
could be as large as one hundred per year, for signals detectable out to hundreds of Mpc 
by the advanced version of LIGO, and because the signal can be accurately predicted using 
general relativity. 

That the signal can be calculated with high accuracy is essential for the measurement 
of the source parameters 0, which include distance, position in the sky, orientation of the 
orbital plane, and the masses and spins of the companions. Loosely speaking, the measured 
signal is passed through a linear filter constructed from the expected signal h(t; 6) and the 
spectral density of the detector noise |J (see below). The signal-to- noise ratio S/N(6) is 
then computed. The expected signal and the signal-to-noise ratio are expressed as functions 
of the vector 6 which collectively represents the source parameters. The actual value of 
these parameters, which we denote 6, is unknown prior to the measurement. When 6 = 
the linear filter becomes the Wiener optimum filter which is well known to yield the largest 
possible signal-to- noise ratio ||. The source parameters can therefore be determined by 
maximizing S/N(0) over a broad collection of expected signals h(t; 6), loosely referred to as 
"templates." 

The gravitational-wave signal can be characterized by a growing amplitude and a phase 
which accumulates nonlinearly with time [lj. The signal undergoes a number iV of oscil- 
lations, varying from 600 to 16 000 depending on the nature of the system (see below), as 
the frequency sweeps through the detector bandwidth. It has been established that it is the 
phasing of the signal which plays the largest role in parameter estimation |7]|],|nj. This is 
because a slight variation in the parameters can quickly cause h(t; 6) to get out of phase 
with respect to the true signal h(t;0), thus seriously reducing S/N(6) from its maximum 
possible value. Therefore a good match between the template's phase and that of the mea- 
sured signal, throughout the N cycles, singles out, to a large extent, the value of the source 
parameters. (Clearly, this is only true for those parameters which affect the phasing of the 
waves, such as the masses and spins of the companions; see below.) 

In principle, the gravitational-wave signal from an inspiralling compact binary can be 
calculated exactly using general relativity (this would require the numerical integration of 
Einstein's equations). In practice, however, one must rely on some approximation scheme. 



It appears appropriate, in this context, to adopt a slow-motion approximation |Tl|], and 
to solve the field equations using a combination of post-Newtonian and post-Minkowskian 
expansions fll2 |. To date, the waveform has been calculated accurately through order (f/c) 4 



(where v is the orbital velocity) beyond the leading-order, quadrupole-formula expression 
fl3HT5|. Leading-order expressions are referred to as "Newtonian;" the waveform is therefore 
known to second-post-Newtonian, or 2PN, order. 
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The detailed expression for the 2PN waveform is complicated: the dependence on the 
various angles (position of the source in the sky, orientation of the detector, orientation of 
the polarization axes) is not simple, and the waves have several frequency components given 
by the harmonics of the orbital frequency (assuming that the orbit is circular [fL6|J17|l ). For 
the purpose of this paper, and following Cutler and Flanagan ]I(J , we shall use a simplified 
expression for the waveform. We shall ignore all post-Newtonian corrections to the wave's 
amplitude, and single out its dominant frequency component at twice the orbital frequency 
IJ. Thus, setting G — c — 1, 

h(t;0) =r- 1 Q(angles)A<(7rA<F) 2/3 cos<l>(t). (1.1) 

Here, r is the distance to the source, Q a function of the various angles mentioned above, 
F(t) the gravitational- wave frequency, and §(t) = J27iF(t)dt the phase. We have also 
introduced the (so-called) chirp mass A4: If /i = mivn^j + m 2 ) is the reduced mass and 
M = mi + m2 the total mass, then 

M = r] 3/5 M, 7] = fi/M. (1.2) 

In Eq. ( |1.1| ) we use the most accurate expression available for the phase function $(£). 



It is determined by the 2PN expression for the frequency sweep |13 
dF 96 
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Apart from the parameters introduced previously, dF/dt also depends on L, the direction 
of orbital angular momentum, and on Si and S%, the spin angular momentum of each 



companion. This dependence is hidden in the "spin-orbit" parameter [18 



ft = ~TK Et 113 K/ M ) 2 + 75 V]L ■ X,, (1.4) 
iZ i=i 



where Xi — Si/ m i i an d the "spin-spin" parameter [[18 



a = ^(-247xi • X2 + 721L ■ X iL ■ * 2 ). (1.5) 

The purpose of this paper is to estimate the anticipated accuracy with which the various 
parameters (such as A4, rj, (3, and a) can be determined during a gravitational- wave mea- 
surement. This analysis differs from that of Cutler and Flanagan []10| in that it incorporates 
terms of 2PN order into the phasing [terms of order (nMF) 4 / 3 in Eq. ( |1.2| )1; their calcu- 
lations were accurate only through 1.5PN order (terms of order 7rMF). Previous analyses 
also include Refs. @gg-gg. 

Our main conclusion is that 1.5PN phasing underestimates the uncertainty in such pa- 
rameters as A4, 7], and (3: 2PN phasing predicts somewhat larger measurement errors. This 
is true even when no attempt is made to determine the spin-spin parameter a. If, however, 
a is also estimated, then the measurement errors become even larger. This is because the 
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TABLE I. Contributions to the accumulated number of wave cycles measured in a 
LIGO/Virgo-type detector. The frequency entering the bandwidth is 10 Hz (seismic limit); the 
frequency leaving is 1000 Hz (system A) (shot noise), and 360 Hz (system B) and 190 Hz (system 
C) (innermost circular orbit). The various contributions correspond to various terms in Eq. (1.3). 
Newtonian: first term within the square brackets; 1PN: second term; tail: Ait^irMF); spin-orbit: 
-P(ttMF); 2PN: (vrMF) 4 / 3 terms with a = 0; and spin-spin: cr(vrMF) 4 / 3 . 



system A B C 

Newtonian 16 050 3 580 600 

1PN 439 212 59 

tail -208 -180 -51 

spin-orbit 17/3 14/3 4/3 

2PN 9 10 4 

spin-spin —2a —3a —a 



number of estimated parameters has increased with respect to the number contained in the 
1.5PN waveform. 



An independent analysis by Cutler and Flanagan [22|] shows that 2PN waveforms are not 
sufficiently accurate for the purpose of parameter estimation: they produce systematic errors 
which are larger than the statistical errors inherent to the measurement process. This is 
because the 2PN waveform fails to remain in phase with the true general-relativistic signal, 



even when the source parameters are exactly matched [[23| -|25ll- To construct templates such 
that the systematic errors will fall below the measurement errors will require an expression 
for the wave's phasing accurate through at least 3PN order [p^-p5|. To achieve such a high 



degree of accuracy is a major challenge for gravitational-wave theorists. 
We shall consider the following three "canonical" binary systems: 

System A: two neutron stars, with mi = m 2 = 1.4M ; 

System B: neutron star and black hole, with mi = 1.4M (the neutron star) and 
m 2 = 1OM (the black hole); and 

System C: two black holes, with mi = m 2 = 1OM . 

For each of these systems Table I summarizes the contribution from each term in Eq. ( |1.3| ) 
to the total number of gravitational- wave cycles received in a LIGO/Virgo-type detector. 

The remainder of the paper is organized as follows. In Sec. II we review the theory of 
parameter estimation, as developed in previous papers by Finn [p6| and Cutler and Flanagan 
||10|| . In Sec. Ill we carry out the calculations for the waveform ( |1.1| ), and compute the 
anticipated uncertainty in the measured values of the source parameters. Our results are 
summarized and discussed in Sec. IV. 



II. PARAMETER ESTIMATION: THEORY 

The theory of detection and measurement of gravitational-wave signals was put on a firm 



statistical foundation, rather similar to that underlying the theory of radar detection [p|,|27 
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This was done by various authors, including Finn p6[ and Cutler and Flanagan |10 |. In this 



section we review the various aspects of the theory which are relevant for our purposes. 

We assume that some criterion has been applied to conclude that a signal originating 
from an inspiralling compact binary has been received by a network of gravitational-wave 
detectors. It is therefore known that a signal of the form h(t; 6) has passed through the de- 
tectors, and we seek to determine the value of the source parameters and the measurement 
error AO = 6 — 0, where 6 denotes the true value. 



Finn |26] has derived an expression for p(6\s), the probability that the gravitational- wave 
signal is characterized by the parameters 0, given that the detector output is s(t) and that 
a signal h(t; 0) — for any value of the parameters — is present. The detector output is 
given by 

s(t) = h(t;0)+n(t), (2.1) 

where n(t) represents the detector noise, assumed to be a stationary, Gaussian random 
process. Finn shows that 

p(0\s) ocp (o) (0)exp[-i(/i(0) -s | h(0)-s)], (2.2) 

where p(°'(0) is the a priori probability that the signal is characterized by 6 (this represents 
our prior information regarding the possible value of the parameters) and where the constant 
of proportionality is independent of 0. 



The inner product (-|-) is defined as follows [10|. The statistical properties of the detector 



noise can be summarized by its autocorrelation function C n (r) = (n(t)n(t + r)), where (•) 
denotes a time average. (It is assumed that the noise has zero mean.) The Fourier transform 
of the autocorrelation function gives the noise spectral density 

/GO 

C n (r)e 2 ^ r dr, (2.3) 
-oo 

which is defined for / > only. The inner product is defined so that the probability for the 
noise n(t) to have a particular realization no(t) is given by p(n = no) oc exp[— (no|no)/2]. It 
is given by 

(, w= ^ ray% , (2 . 4 ) 

where g(f) is the Fourier transform of g(t), 

9(f) = / 9(t)e 2mft dt; (2.5) 



an asterisk denotes complex conjugation. 

We define p, the signal-to-noise ratio associated with the measurement, to be the norm 
of the signal hit] 6), 

P 2 = (h\h) = 4 f l 4r£df, (2.6) 



o 



Sn(f) 
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evaluated at 6 = 6, the true value of the source parameters. This is the largest possible 
value of the signal-to-noise ratio, for h*(f;0)/ S n (f) is just the Fourier transform of the 
Wiener optimum filter ||. 

In a given measurement, characterized by the particular detector output s(t), the true 
value of the source parameters can be estimated by locating the value 6 at which the prob- 
ability distribution function ( |2.2p is a maximum. This is the so-called maximum-likelihood 
estimator ||. In the limit of large signal-to-noise ratio, to which we henceforth specialize, 
p{0\s) will be strongly peaked about that value. We now derive a simplified expression for 
p{0\s) appropriate for this limiting case. 

We first assume that p(°\0) is nearly uniform near 6 = 0. This indicates that the prior 
information is practically irrelevant to the determination of the source parameters; we shall 
relax this assumption below. Then, denoting = (h(0) — s\h(0) — s), we have that £ is 
minimum at 6 = 0. It follows that this can be expanded as 

ae) = m + ^a b (e)Ae a Ae b + --- J (2.7) 

where ", a" denotes partial differentiation with respect to the parameter 8 a , and A9 a = 
9 a — 9 a ; summation over repeated indices is understood. We assume that p is sufficiently large 
that the higher-order terms can be neglected. Calculation yields ^ a b = (h ta b\h — s) + (h ya \h t b), 
and we assume once more that p is large enough that the first term can be neglected (see 
Cutler and Flanagan [plj for details). We arrive at 



p(9\s) ocp (o) (0)exp[-ir afe A0 a A(9 6 ], (2.8) 

where 

T ab = (h a | h b ), (2.9) 



evaluated at = 6, is the Fisher information matrix |27 |. 

We therefore see that in the limit of large signal-to-noise ratio, p{0\s) takes a Gaussian 
form. From Eq. ( [2.8|) it can be established that the variance-covariance matrix E a6 is given 
by 



Here, (■) denotes an average over the probability distribution function (|2.8|) , and T 1 is the 
inverse of the Fisher matrix. We define the measurement error in the parameter 9 a to be 



a a = ((A6 a ) 2 } 1/2 = Vf^ (2.11) 

(no summation over repeated indices), and the correlation coefficient between parameters 
9 a and 9 b as 

by definition each c ab must lie in the range (—1,1). 
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Cutler and Flanagan [10 



Eq. 



have shown that in the limit of large signal-to-noise ratio, 
is valid even when p(°'(6) is not uniform near = 0. In such cases the prior 



information plays an important role in the determination of the source parameters. The ex- 
ponential factor is still peaked at = 0, but no longer represents the maximum-likelihood 
estimate, and the full probability distribution function p(0\s) may not be a Gaussian. 

For simplicity, and following Cutler and Flanagan ||10|| , we shall restrict attention to cases 
such that (0) is a Gaussian, given by 



p (o) (0) ocexp 



i r (o), 

~2 L oh \ 



a ){6 b -6 h 



(2.13) 



Then p(0\s) will also take a Gaussian form, and the new variance-covariance matrix will be 
given by 



(r + r<°>) 



-i 



(2.14) 



It should be noted that in general, p{0\s) will be peaked at a value (0) which differs both 
from and 0. 



III. PARAMETER ESTIMATION: CALCULATIONS 



We proceed with the calculation of the Fisher information matrix, Eq. ( |2.9[ ), for 
gravitational- wave signals of the form (|1 . 1|) , and for gravitational- wave detectors of the 
LIGO/ Virgo type. For such detectors the anticipated noise spectral density can be approx- 
imated by the analytic expression [|1(| 



^(/) = ^ [(/o//) 4 + 2 + 2(/// ) 



(3.1) 



where Sq is a normalization constant irrelevant for our purposes, and /o the frequency 
at which S n (f) is minimum; we set Jq = 70 Hz, which is appropriate for advanced LIGO 
sensitivity ||. To mimic seismic noise we assume that Eq. ( |3.1| ) is valid for / > 10Hz only, 
and that S n (f) = oo for / < 10Hz. 

First, we integrate Eq. (|1.3| ) to obtain expressions for <&(F) and t(F), respectively the 
phase and time as functions of gravitational-wave frequency. (Throughout this section we 
shall distinguish between F, the function of time describing the frequency sweep, and /, the 
Fourier-transform variable.) Expanding in powers of (nMF) 1 ^ 3 and truncating all expres- 
sions to 2PN order, we obtain 



d>(F) 



- y^MF) 



-5/3 
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1 + 1 n 
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n h 7T 
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a 



{■kMF) 
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(3.2) 



where C is (formally) the value of $ at F = oo, and 



t(F) = t c 



256 



M(nMF) 



-8/3 



/ 3058673 5429 

+ 2 1 

V 1016064 1008 



4/743 

1 + - 

3 V336 

617 



—rj) (nMF) 2/3 - - (4tt - (3) (rtMF) 
4 / 5 



(3.3) 
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where (formally) t c = t(oo). Of course, the signal cannot be allowed to reach arbitrarily 
high frequencies; it must be cut off at a frequency F = F; L corresponding to the end of the 
inspiral. We put ixMFi = (M/r^) 3 / 2 = 6~ 3 / 2 ; = 6M is the Schwarzschild radius of the 
innermost circular orbit for a test mass moving in the gravitational field of a mass M |28[ . 
Next, we take the Fourier transform of Eq. (EI) and calculate h(f) = J hfye^^dt. It 



is sufficient to estimate h(f) using the stationary phase approximation, according to which 



l{t)e im dt 



2vu 



1/2 



g(to)e 



i(j>(to) 



(3.4) 



if g(t) varies slowly near t = t where the phase has a stationary point: 0'(^o) = (a 
prime denotes differentiation with respect to t). Substituting Eqs. (|1.1|) and ( |3.2|) into (p.4|), 
discarding the irrelevant negative-frequency component, and neglecting all post-Newtonian 
corrections to the amplitude of h(f), we obtain 



h(f) = Af- 7/6 e^ f \ 
where A oc M. 5 > 6 Q( angles )/r, and 



(3.5) 



v>(/) = 2nft c - 0c - ^ + ^(^fr 5/3 



20 /743 11 > , r . ., 
1 + — I — + — v) (ttM/) 2 / 3 
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1016064 1008 



144 



^-^(ttM/) 4 / 3 



(3.6) 



We have introduced the parameter e = 1. This gives us the freedom, for future use, of 
removing the 2PN terms from ip(f) by setting e = 0. 

We now substitute Eq. ( |3.5| ) into Q2.6| ) and calculate the signal-to-noise ratio. We readily 
obtain 

p 2 = 20^ 2 5 -Vo~ 4/3 /(7), (3.7) 
where the integrals I(q) represent various moments of the noise spectral density: 



Hi) 



— 7;dx, 
1/7 x- 4 + 2 + 2x 2 



(3.8) 



where Xi = fi/fo = (6 3 ^ 2 vrM/o) _1 is the frequency cutoff. 

As the next step toward the computation of the Fisher matrix, we calculate the deriva- 
tives of h(f) with respect to the seven parameters 



6 = (In A, f t c , 4> c , \nM, In 77, (3, a) 



(3.9) 



We obtain 



hi = h, 

~h 2 = 2m{f/f )h, 



h 



,3 



—ih, 
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57 

h A = - — (7rMf)- 5/3 (l + A,v 2 - B 4 v 3 + C,v 4 )h, (3.10) 

h, 5 = -l-^MfY^iAsv 2 - B 5 v 3 + C 5 v 4 )h, 
96 



64 

where v = (ttM/) 1 / 3 . We also have defined 

4 /743 11 
4 ~ 3 V336 + ~J V 



S 4 = -(47r-/5), (3.11) 

„ /3058673 5429 617 2 

C 4 = 2e 1 77 H r/ 2 - a 
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and 



743 33 

Ak = 77, 

5 168 4 
27 

fi 5 = -(4vr-/3), (3.12) 

10 /3058673 5429 617 2 

C 5 = 18e 7? tT - a 

V 1016064 4032 ' 96 ' 

Finally, the components of T can be obtained by evaluating the inner products (h ta \h t b) 
using Eq. ( |2.4|) . The r a b's can all be expressed in terms of the parameters 6, the signal-to- 
noise ratio p, and the integrals I(q). The expressions are too numerous and lengthy to be 
displayed here. As illustrating examples, we quote 

r ia = 5 laP \ (3.13) 

and 

r 46 = - J^ 3/ VWo)- 7/; V(i4) + A 4 J(i2) 

-y3 4 J(ll) + C 4 J(10)]p 2 , (3.14) 

where J(q) = I(q)/I(7). We note that even though the h a 's are expressed as truncated 
post-Newtonian expansions in Eq. ( |3.10|) , they must be treated as exact when computing T. 
This is to ensure that the eigenvalues of the Fisher matrix are always positive definite. 
The variance-covariance matrix S afe can now be obtained from Eq. ( 2.14|) , and the mea- 



surement errors and correlation coefficients computed from Eqs. ( [2.11 ) and (|2.12| ). Before 



doing so, however, we must first state our assumptions regarding the prior information 
available on the source parameters. 

We take advantage of the fact that the dimensionless spin parameters, \i an d X2, must 
necessarily be smaller than unity. (This upper bound is strict for black holes, but only 
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TABLE II. Measurement errors and correlation coefficients for the estimation of six parameters 
(<t is not estimated), assuming p = 10 and (3 = 0. The first column indicates whether or not prior 
information was included in the calculation. The second column gives the value of e, introduced 
in Eq. (3.6); e = 1 represents 2PN phasing. Then follows, in more suggestive notation, 02//0 (in 
msec), <7 3 (in radians), 04, 05, ae, c 45 , c 46 , and c 56 (all dimensionless) . 



prior 


e 


At c 


A(f> c 


AM/M 


Arj/r] 


A/3 


c Mr > 


c Mf3 




System 


A (two 


neutron stars): 














yes 


1 


1.07 


2.94 


0.036 % 


0.279 


1.33 


-0.989 


0.994 


-0.999 


no 


1 


1.08 


2.97 


0.037 % 


0.282 


1.35 


-0.989 


0.994 


-0.999 


no 





1.13 


4.09 


0.034 % 


0.243 


1.24 


-0.988 


0.993 


-0.999 


no 


-1 


1.16 


4.96 


0.032 % 


0.213 


1.15 


-0.986 


0.992 


-0.999 


System 


B (neutron star and black hole): 












ye& 


1 


1.72 


2.27 


0.218 % 


503 


2.29 


-0 993 


996 


-0 999 


no 


1 


1.76 


2.32 


0.226 % 


0.523 


2.38 


-0.993 


0.996 


-0.999 


no 





2.04 


6.24 


0.191 % 


0.386 


1.99 


-0.990 


0.994 


-0.999 


no 


-1 


2.20 


8.68 


0.171 % 


0.306 


1.76 


-0.988 


0.993 


-0.999 


System 


C (two 


black holes): 














yes 


1 


1.50 


2.19 


0.54 % 


1.46 


8.19 


-0.946 


0.956 


-0.999 


no 


1 


2.40 


4.99 


1.96 % 


5.50 


30.8 


-0.996 


0.997 


-0.999 


no 





3.53 


9.27 


1.42 % 


3.16 


19.5 


-0.992 


0.994 


-0.999 


no 


-1 


4.01 


14.7 


1.21 % 


2.22 


14.9 


-0.989 


0.992 


-0.999 



approximate for neutron stars.) It follows from Eqs. ( |1.4| ) and (|1.5|) that must be smaller 
than approximately 8.5, and that |er| must be smaller than approximately 5.0. Following 
Cutler and Flanagan flOfl , we crudely incorporate this information into our calculations by 
taking 

p<°>(0) oc exp[-i(/3/8.5) 2 - |(a/5.0) 2 ]. (3.15) 
We consider all other parameters to be unconstrained |30| . 



IV. RESULTS AND DISCUSSION 



Equation (|3.12|) implies that the Fisher matrix is block diagonal. The parameter 9 l 



hi A is therefore entirely uncorrelated with the other parameters, and we find o"i = AA/A = 
1/p, c la = 0, in all cases. We shall no longer be concerned with this parameter. 

The results concerning the other parameters are displayed in Tables II and III. All 
calculations were carried out assuming p = 10, and that the companions are spinless, so 
that f3 = a = 0. 

To obtain the results of Table II we have estimated only six (including A) of the seven 
parameters, leaving a out. In effect, we have truncated the original Fisher matrix to a 
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smaller, 6x6, matrix. This amounts to assuming before measurement that the spin-spin 
parameter must be very small; equivalently, this assumption can be implemented by taking 
p (0) O) to be very strongly peaked at a = 0. 

For each of the three canonical systems, the first line of Table II displays the measure- 
ment errors and correlation coefficients as calculated using 2PN phasing (e = 1) and the 
prior probability distribution function ( ft. 15 ). The second line shows the same quantities 
calculated without utilizing the prior information. We notice that the prior information 
makes virtually no difference for systems A and B, but is very significant for system C. 

The third line of Table II displays the measurement errors and correlation coefficients 
assuming no prior information and 1.5PN phasing (e = 0). Our values agree with those of 
Cutler and Flanagan [|l0]|3l|. We notice that the errors calculated using 1.5PN phasing are 



always larger hi t c and <p c , and smaller for M., T], and (3, than those calculated using 2PN 
phasing. Thus, the measurement errors on the masses and spins are underestimated when 
evaluated using the less accurate 1.5PN phasing. 

This can be explained with a simple argument. In Eq. ( |3.6| ), the 1PN and 2PN terms [of 
order (nMf) 2 ^ 3 and (7rM/) 4 / 3 respectively] combine, when e = 1, so as to reduce the relative 
importance of the 7rMf term, when compared to the situation when e = 0. In other words, 
the relative contribution to the total number of wave cycles coming from the nMf term is 
less for 2PN phasing than it is for 1.5PN phasing (see Table I). It is therefore expected that 
2PN phasing will produce larger measurement errors for f3, since all information about (3 
comes from the nMf term. But because (3 is strongly correlated with both M. and rj, it 
follows that these parameters will also come with larger measurement errors. This is indeed 
what is observed. It is amusing to test this explanation by artificially setting e = — 1 in 
our calculations, which we do in the fourth line of Table II. The argument suggests that the 
errors in M., rj, and (3 should all decrease with respect to the values calculated using 1.5PN 
phasing, since the relative importance of the nMf term is now increased. This is indeed 
what the results show. 

To obtain the results of Table III we have estimated all seven parameters, including 
both A and a, and used 2PN phasing. For each of the three systems, the first line of 
Table III displays the measurement errors and correlation coefficients calculated using the 
prior probability distribution function ( ft.!5| ). We notice that the measurement errors are 
all significantly larger than those displayed in Table II; this is expected from the fact that 
we are now estimating a larger number of parameters. 

It is interesting to ask how the measurement errors increase as the number of estimated 



parameters increases. In Ref. ||10|| , Cutler and Flanagan initially estimate only five of their 
six parameters, leaving (3 out. When they next include f3 in their calculations, they find 
that the measurement errors on M. and fi increase by a factor of order 10. In this paper, on 
the other hand, we have initially estimated only six of our seven parameters, leaving a out. 
When we next include a in our calculations, we find that the measurement errors on M, //, 
and j3 only increase by a factor of order unity. Thus, the inclusion of a in the calculation has 
less dramatic consequences than the inclusion of (3. This confirms a conjecture formulated 
by Cutler and Flanagan [10 at the end of their Sec. III. That this is so is largely due to the 
importance of prior information in the estimation of a. 

In the second line of Table III we display the results obtained when the prior information 
is not included in the calculations. We notice that for all systems, the prior information 
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TABLE III. Measurement errors and correlation coefficients for the estimation of all seven 
parameters, assuming p = 10, (3 = 0, and a = 0. The first column indicates whether or not prior 
information was included in the calculation. The notation is similar to that of Table II, and e = 1 

in all cases. 

prior At c Act> c AM/M A 77/77 A/3 Act c m ^ c M/3 c Ma cP a 
System A (two neutron stars): 

yes 1.28 13.3 0.047 % 0.507 1.77 4.79 -0.956 0.996 -0.648 -0.964 0.835 -0.660 
no 2.72 46.9 0.120 % 1.578 4.53 17.3 -0.991 0.999 -0.952 -0.993 0.984 -0.955 
System B (neutron star and black hole): 

yes 2.54 23.6 0.280 % 0.873 3.02 4.74 -0.959 0.997 -0.630 -0.969 0.817 -0.650 
no 7.52 95.9 0.813 % 3.10 8.98 19.4 -0.993 0.999 -0.961 0.995 0.986 -0.964 
System C (two black holes): 

yes 2.22 10.4 0.55 % 1.51 8.22 4.81 -0.849 0.920 0.191 -0.984 0.257 -0.081 
no 17.0 179 7.23 % 30.7 149 74.6 -0.995 0.998 -0.962 -0.999 0.984 -0.978 



indeed plays a very important role. In fact, we see that Eq. (|3.15|) provides nearly all of 
the information regarding the spin-spin parameter a. This explains why the measurement 
error on a is always nearly equal to 5.0, and its correlation coefficient with other parameters 
significantly smaller than unity. Of course, these results only apply to gravitational-wave 
measurements with p = 10. To bring the error on a well below the a priori constraint 
cr < |5.0|, say Ao < 3, the measurement would require a signal-to-noise ratio larger than 
approximately 45 for system A, 50 for system B, and 110 for system C. 

We conclude with the following remark. It is clear that the results displayed in Tables II 
and III depend on a fairly large number of simplifying assumptions, and that a more careful 
treatment might produce somewhat different numbers. These assumptions include: (i) the 
simplified form ( |1.1| ) for the waveform; (ii) the neglect of (not yet calculated) higher-order 
terms in the post-Newtonian expansion ( |1.3|) ; (iii) the neglect of 0(lp) corrections in the 
expression (|2.8j), ( |2.9| ) for p(6\s); (iv) the analytic model ( |3.1| ) for the noise spectral density; 
and (v) our rather crude incorporation of the prior information. We shall leave for future 
work the difficult task of carefully examining the effect of these assumptions on our results. 
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